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摘 要 : 利用 陕西 省 某 一 风电 场 区 域内 的 观测 资料 ,分 析 了 该 风电 场 的 风速 规律 ,并 引入 最 优 训 练 
期 方案 ,研究 利用 线性 回归 方法 建立 风电 功率 预报 模型 的 可 行 性 。 结 果 表 明 : 该 风电 场 区 域 ,不 同 
高 度 的 风速 及 其 高 度 间 的 风速 差异 均 表现 出 最 大 值 出 现在 夜间 ,最 小 值 出 现在 白天 ,从 低层 到 高 
层 的 风速 日 变化 趋势 一 致 的 特征 。 一 日 中 ,风速 与 风电 功率 在 09:00 ~ 17:00 时 段 的 相关 系数 明显 
小 于 其 它 时 段 。 按 照 风 速 是 否 大 于 5m's "将 训练 期 观测 样本 分 为 2 组 ,可 以 明显 改善 风速 与 风电 
功率 的 回归 关系 。 以 风机 轮 载 高 度 处 的 风速 作为 预报 因子 ,并 引入 风电 功率 与 风速 之 闯 相关 系数 
的 日 变化 规律 、 以 及 不 同 风 速 量 级 下 风速 与 风电 功率 之 间 回 归 关 系 的 差异 性 ,采用 最 优 训练 期 方 
案 和 一 元 线性 回归 方法 建立 的 风电 功率 预报 方程 ,具有 预报 误差 小 和 最 优 训练 期 短 的 特点 ,满足 


实际 业务 需求 。 
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随 着 全 球 气候 变 暖 和 能 源 紧张 等 问题 的 日 益 
严峻 ,针对 风能 等 可 再 生 能 源 的 研究 越 来 越 受 到 重 
视 。 大 容量 的 风电 并 入 电网 会 影响 电力 系统 的 稳 
定 运行 ,对 风电 合理 调配 利用 、 风 电 市 场 管理 及 商 
业 运营 等 带 来 严峻 挑战 于 ,因此 对 风电 功率 的 预测 
研究 提出 了 迫切 需求 。 

欧洲 一 些 国家 的 风能 预报 工作 起 步 较 早 ,目前 
已 经 有 不 少 比较 成 熟 的 风能 预报 系统 ,如 丹麦 的 
PREDIKTOR 和 WPPT 美国 的 EWIND 加 拿 大 的 
WEST .德国 的 PREVIENTO .西班牙 的 LOCALPRED 
和 REGIOPRED 等 。 我 国 风电 场 风能 预报 工作 开 
展 较 晚 ,科研 单位 主要 有 中 国电 力 科 学 研究 院 .中 
国 气 象 局 国家 气候 中 心 .华北 电力 大 学 和 清华 大 学 
等 由 。 风 电功率 预测 的 方法 有 两 种 ,一 种 是 利用 实 
况 资料 ,通过 外 推 的 方法 进行 预报 ,如 自 回归 滑动 
平均 法 5 .时间 序 列 法 5、 卡尔 曼 滤 波 法 以 及 人 工 
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神经 网 络 方法 等 ,但 是 这 种 方法 计算 结果 的 准确 
率 和 时 效 性 较 低 。 男 一 种 是 依靠 数值 天 气 预 报 制 
作风 电功率 预报 ”中 ,预测 时 间 可 以 达到 24 ~72h 
或 者 更 长 ,但 是 预报 精度 会 随 着 预报 时 效 的 增加 而 
降低 。 为 了 提高 预测 精度 ,也 有 研究 者 提出 了 组 合 
预测 方法 ,该 方法 减少 了 较 大 误差 点 ,但 是 模型 复 
Ze RABE. 

在 众多 的 统计 分 析 方 法 中 ,线性 回归 方法 是 最 
第 用 的 方法 之 一 ,具有 所 需 资料 少 、 计 算 过 程 简便 
的 特点 。 利 用 气温 预报 值 与 观测 值 之 间 的 一 元 线 
性 回归 关系 建立 的 误差 订正 方程 ,可 以 改善 模式 的 
气温 预报 质量 ,并 被 气象 部 门 广泛 应 用 于 气温 预报 
业务 中。 基于 一 元 线性 回归 方法 , 匡 启 树 等 ' 引 设 
计 了 最 优 训练 期 方案 对 ECMWF 模 式 的 日 最 高 ( 低 ) 
气温 进行 预报 ,显著 提高 了 气温 预报 质量 。 那 么 ， 
利用 风速 与 风电 功率 输出 之 间 的 统计 关系 ,通过 线 
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性 回归 方法 建立 的 风电 功率 预报 模型 是 否 也 具有 
业务 应 用 价值 ? 另外 ,风速 是 影响 风能 资源 开发 利 
用 的 重要 因子 ,研究 不 同 区 域 风 速 的 变化 规律 ,对 
提高 风电 场 的 功率 预测 水 平和 做 好 风能 资源 评估 
具有 重要 意义 “”。 本 文 将 利用 风电 场 的 观测 资 
料 ,分 析 风 电场 区 域内 的 风速 变化 规律 ,并 引入 最 
优 训练 期 方案 ,研究 利用 线性 回归 方法 建立 风电 功 
率 预 报 模型 的 可 行 性 。 


1 资料 和 处 理 


研究 资料 来 自 2016 年 1 月 1 日 ~2018 年 12 月 
31 日 陕西 省 某 一 风电 场 区 域内 1 台 测 风 塔 和 33 A 
风机 的 逐 5 min 观 测 ,由 于 资料 保密 原因 ,省 略 该 风 
电场 的 具体 地 理 位 置信 息 。 测 风 塔 的 观测 要 素 包 
括 风速 .风向 .温度 和 气压 ,其 中 ,风速 和 风向 的 测 
量 高 度 分 别 是 80 m、50 m,30 m 和 10 m, 温 度 和 气压 
的 测量 高 度 均 为 10 m。 风 机 的 观测 要 素 包括 风速 
和 有 功 功率 ,风机 轮 载 高 度 为 80 ms 

由 于 天 气 或 机 器 故障 等 原因 , 测 风 塔 和 风机 的 
观测 数据 会 有 不 合理 值 出 现 , 为 了 保证 研究 结论 的 
可 靠 性 ,有 必要 剔除 异常 值 汪 。 对 测 风 塔 观测 资料 
的 质量 控制 包括 以 下 方面 :(1) 风 速 和 风向 的 观测 
值 范 围 分 别 在 0~40 m- s~ FAO. ~360° 之 间 。(2) 风 
WEAF 5 ms"!' 且 风速 和 风向 连续 6 无 变化 的 情况 
视 为 缺 测 ,风速 小 于 0.4 ms"' 且 连续 3 以 上 无 变化 
的 情况 视 为 缺 测 。(3) 相 隔 高 度 在 1~20m 和 21~ 
40m 条 件 下 的 风速 差 值 分 别 小 于 $ mes #10 m.s 
(4) 当 风速 大 于 $m:s 时 ,风速 的 标准 偏差 小 于 
10 ms'。(5) 高 度 风向 差 小 于 22.5°。(6) 小 时 平均 风 
速 ( 温 度 ) 变 化 小 于 10 ms'(5 C),3 hE) F 1 
hPa。 对 风机 观测 资料 的 质量 控制 包括 以 下 方面 : 
(1) 风 速 大 于 5 ms 而 有 功 功率 小 于 50kW 的 记录 
值 子 以 剔除 。(2) 风 速 和 有 功 功率 分 别 大 于 切 出 风速 
20 ms" 和 额定 功率 1500 kW 的 记录 值 耶 以 剔除 。 


2 风电 场 风速 规律 研究 


利用 测 风 塔 和 风机 的 逐 5 min 观测 资料 ,统计 得 
到 风机 有 功 功 率 与 风机 风速 测 风 塔 的 80 m 50 m 
30 m Fil 10 m 风 速 的 相关 系数 依次 为 0.88 , 0.84, 
0.83 .0.81 和 0.75( 均 通过 了 a = 0.01 的 显著 性 水 平 检 
验 ), 有 明显 的 正 相 关 关 系 , 与 10m 气 温和 10m 和 气压 


的 相关 系数 分 别 为 -0.12 和 -0.06 ,相关 性 不 明显 ， 
可 见 , 风 速 是 影响 风电 功率 的 最 主要 因子 。 为 了 了 
解 风 电场 所 在 区 域内 的 风速 变化 和 风 廊 线 特征 ,下 
面 对 测 风 塔 不 同 高 度 的 风速 观测 结果 进行 分 析 。 

图 1 为 2016 年 1 月 1 日 ~2018 年 12 月 31 日 测 
风 塔 的 80 m、50 m、30 m 和 10 m 不 同 高 度 各 季节 平 
均 风速 的 日 变化 。 从 图 1 来 看 ,春季 秋季 和 冬季 的 
风速 日 变化 特征 相近 ,各 高 度 的 平均 风速 从 凌晨 开 
始 逐 渐 减 少 ,在 09:00 ~ 10:00 左 右 达 到 极 小 值 ,然后 
开始 增 大 ,到 14:00 左 右 达 到 极 大 值 ,然后 再 次 逐渐 
减 小 ,在 19:00 左 右 再 次 达到 极 小 值 ,之 后 又 逐渐 增 
大 ,直至 凌晨 23:00~ 次 日 01:00 左 右 达 到 极 大 值 。 
夏季 ,80m.50m 和 30m 各 高 度 的 平均 风速 从 凌晨 
开始 逐渐 减少 ,在 08:00 ~ 09:00 左 右 达 到 极 小 值 , 然 
后 开始 逐渐 增 大 ,直至 凌晨 02:00 左 右 达 到 极 大 值 ， 
而 10 mm 平均 风速 的 日 变化 特征 与 其 它 季节 相似 。 
整体 上 讲 ,夜间 的 风速 大 于 白天 ,风速 极 小 值 出 现 
在 上 午 和 傍晚 , 极 大 值 出 现在 凌晨 和 下 午 ,页 莹 莹 
等 研究 了 我 国 风蚀 区 的 风速 日 变化 特征 ,发现 部 
分 山区 的 站 点 日 最 大 风速 出 现在 夜间 ,与 本 文 的 研 
究 结论 一 致 。 另 外 ,风机 有 功 功 率 的 最 小 值 出 现在 
上 午 ,最 大 值 出 现在 凌晨 ,与 图 1 给 出 的 测 风 塔 观 
测 到 的 风速 的 日 变化 特征 有 一 些 差异 (如 ,在 下 午 
14:00 左 右 没有 出 现 有 功 功 率 的 极 大 值 ) ,但 是 与 风 
机 站 观测 到 的 风速 的 日 变化 特征 一 致 (图 略 )。 

风速 的 垂直 分 布 会 影响 风电 功率 输出 ,如 果 能 
将 风速 廊 线 随时 间 的 变化 规律 引入 风速 和 风电 功 
率 预报 模型 ,将 有 利于 减少 预报 误差 "。 图 2 给 出 
了 各 季节 不 同 高 度 间 风速 差异 的 日 变化 曲线 。 从 
图 2 来 看 ,各 季节 不 同 高 度 间 的 风速 差异 均 表 现 出 
由 夜间 到 白天 逐渐 减 小 ,由 白天 到 夜间 逐渐 增 大 的 
趋势 ,最 小 值 出 现在 上 午 10:00 ~ 11:00 左 右 。 白 天 
的 不 同 高 度 间 的 风速 差异 小 于 夜间 ,这 与 白天 大 气 
层 结 处 于 不 稳定 或 中 性 状态 ,大 气 汕 流 混合 较 强 ， 
有 利于 减 小 各 层 间 的 风速 差异 ,而 夜间 大 气 层 结 稳 
定 , 沸 流 微 纶 ,动量 不 宜 下 传 ,使 得 各 层 间 的 风速 差 
异 较 大 有 关 ””。80 m 与 10 m 高 度 间 的 风速 差异 全 
年 始终 保持 最 大 ,50 m 与 10 m 高 度 间 的 风速 差异 次 
之 。 但 是 不 同 高 度 间 的 风速 差异 与 它们 之 间 的 高 
度 差 并 不 完全 成 正比 ,如 ,30 m 与 10 m 高 度 间 的 风 
速 差异 大 于 与 80 m 高 度 间 的 风速 差异 , 即 , 近 地 面 
到 30 m 高 度 的 风速 减 小 幅度 大 于 30 m 到 80 mm 高度 
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分 析 及 风电 功率 预报 方法 研究 


一 一 10 m 风 速 
一 旦 一 30 m 风 速 
—S— 50 m 风 速 


风速 /ms-! 


1 测 风 塔 各 季节 不 同 高 度 平均 风速 的 日 变化 


Fig.l Diurnal variation of seasonal average wind speed at different heights of wind tower 


的 风速 减 小 幅度 ,这 与 该 风电 场 区 域内 下 垫 面条 件 
对 风速 垂直 变化 的 影响 较 大 有 关 ? 。 

文献 [13] 研 究 的 张 北 和 吉林 两 个 地 区 风电 场 
的 风速 廓 线 分 布 形式 与 本 文 研 究 的 陕西 地 区 风电 
场 相似 ,但 是 各 高 度 的 风速 日 变化 特征 与 本 文 的 研 
究 结 果 有 一 些 差异 ,这 与 各 风电 场 的 下 垫 面 条 件 和 
低层 大 气 层 结 状态 有 一 定 差 异 有 关 。 


3 ”风电 功率 预报 


3.1 方法 

由 于 风电 功率 输出 与 风速 的 相关 性 较 好 ,与 温 
度 和 和 气压 场 的 相关 性 较 弱 ,并 且 孙 川 永 呈 的 研究 也 
已 证 明 温度 场 和 气压 场 的 引入 并 不 能 改进 风电 功 
率 的 预报 结果 ,风向 的 引入 理论 上 虽然 对 风电 功率 
预报 有 重要 意义 ,但 是 预报 结果 并 不 理想 。 本 节 将 
以 风速 作为 预报 因子 ,研究 利用 线性 回归 方法 建立 
风电 功率 预报 模型 的 可 行 性 ,以 及 如 何 选 取 预 报 因 
子 和 设计 训练 期 方案 ,预报 时 效 为 24b 内 逐 1h 预 
报 , 起 报时 间 为 08:00( 北 京 时 间 )。 
3.1.1 线性 回归 方法 在 一 个 给 定 的 风机 站 上 ,对 
于 某 一 预报 时 效 的 有 功 功 率 , 建 立方 程 : 


C,= da,F., +b, (1) 
i=l 


式 中 : C, 为 预报 值 , 已, 为 第 ;个 预报 因子 ,o, WB i 
个 预报 因子 的 回归 系数 权重 ,b 为 常数 项 ,mm 为 参与 集 
成 的 因子 总 数 ,1 为 预报 时 效 ,t=1, 2, …, 24h, 式 中 a, 
和 65 通过 最 小 二 乘法 原理 计算 得 到 。 

3.1.2 最 优 训练 期 方案 采用 滑动 训练 期 ,训练 样 
本 为 预报 日 之 前 NN 日 的 历史 资料 ,分 别 取 ,N=4， 
5, …, 365 d, 逐 日 滑动 建立 风电 功率 预报 模型 ,再 假 
设 预报 日 的 风速 观测 值 为 预报 值 , 代 入 该 预报 模型 
制作 风电 功率 预报 ,最 后 计算 2017 年 1 月 1 日 ~ 
2018 年 12 月 31 日 风电 功率 预报 的 平均 绝对 误差 
(MAE) ,平均 绝对 误差 最 小 值 对 应 的 NN 即 为 最 优 训 
练 期 。MAE 的 计算 公式 为 : 


ix 
== xl, 2 
MAE nd xil (2) 


KEP yA xP AS i RR HR AU EL, n 
为 参与 检验 的 总 预报 次 数 。 

3.1.3 建 模 方案 I 利用 风机 的 风速 和 有 功 功 率 观测 
资料 ,基于 一 元 线性 回归 方法 和 最 优 训练 期 方案 , 设 
计 了 方案 1.1 ~ 1.4( 表 1) , 先 对 每 台风 机 站 进行 功率 
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tió wA 
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时 间 
图 2 风电 场 区 域内 风 廊 线 的 季节 平均 日 变化 特征 


Au80~50 为 80 m 与 50 m 高 度 间 的 平均 风速 差异 ,Au80~30 为 80 m 与 30 m 高 度 间 的 平均 风速 差异 ， 
依 此 类 推 ,Au80~10、Au50~30、Au50~10、Au30~10 分 别 为 相应 高 度 间 的 平均 风速 差异 


Fig.2 Diurnal variation of seasonal average wind profile in wind farm area 


Au80-u50 is the average difference of wind speed between height of 80 meters and 50 meters ,Au80-u30 is the average difference of 


wind speed between height of 80 meters and 30 meters, in the same way, Au80-u10, Au50-u30, Au50-u10 and Au30-ul0 are the 


average difference of wind speed between corresponding heights 


预报 ,然后 累加 后 除 以 风机 台数 得 到 整个 风电 场 的 
平均 功率 ,其 中 ,方案 1.3 和 1.4 按 照 风 速 是 否 大 于 5 
ms ' 将 训练 期 样本 分 为 2 组 后 分 别 建立 预报 方程 。 
当 风 电功率 预报 值 <0 kW 或 者 >1 500 kW 时 ,分 别 今 
风电 功率 预报 值 为 0kW 或 者 1 500 kW. 


表 1 建 模 方案 I 的 方案 设计 


Tab.1 First scheme designs of wind power prediction 


方案 设计 方案 1.1 方案 1.2 方案 1.3 方案 1.4 
是 否 分 预报 时 次 建 模 是 T E 5 
是 否 按照 风速 量 级 对 训练 a= 2 
期 样本 进行 分 组 后 建 模 时 = = 


3.1.4 ZRI 增加 了 其 他 高 度 风 速 作 为 预报 
因子 ,采取 多 元 线性 回归 方法 和 最 优 训练 期 方案 建 
立 风 电功率 预报 模型 。 虽然 风电 功率 与 风机 站 风 
速 的 相关 系数 必然 大 于 与 测 风 塔 不 同 高 度 处 风速 


的 相关 系数 ,但 是 风机 站 与 测 风 塔 的 风速 观测 不 在 
同一 地 点 ,因此 , 建 模 方案 开 没 有 将 风机 站 的 风速 
作为 预报 因子 ,而 是 利用 33 个 风机 站 的 平均 功率 和 
测 风 塔 的 80 m、50 m、30 m、10 m 高 度 风速 的 观测 资 
料 , 将 整个 风电 场 的 平均 功率 作为 一 个 整体 进行 预 
报 。 在 方案 1.1 的 基础 上 , 建 模 方案 开 根 据 预报 因子 
的 不 同 ,设计 了 方案 2.1~2.4( 表 2)。 同样 , 当 风 电 功 
率 预报 值 < 0 kW 或 者 > 1 500 kW 时 ,分 别 令 风 电 功 


表 2 建 模 方 案 II 的 方案 设计 


Tab.2 Second scheme designs of wind power prediction 


方案 设计 80m 50m 30 m 10 m 
预报 因子 风速 风速 风速 风速 
方案 2.1 Vv 

方案 2.2 M v 

方案 2.3 Vv Vv Vv 

方案 24 
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率 预 报 值 为 0kW 或 者 1500kW。 


3.2 结果 分 析 

3.2.1 建 模 方案 I 的 检验 ”对 于 线性 回归 方程 ,样本 
量 的 大 小 直接 影响 着 线性 回归 系数 的 计算 结果 ,但 
是 样本 量 既 不 是 越 大 越 好 ,也 不 是 越 小 越 好 。 最 优 
训练 期 日 数 的 研究 是 为 了 在 预报 日 之 前 的 历史 观 
测 资 料 中 选取 最 合适 的 样本 ,以 建立 风电 功率 与 风 
速 的 最 优 线性 回归 方程 ,从 而 最 大 限度 地 降低 风电 
功率 的 预报 误差 。 从 方案 1.1 ~ 1.4 的 最 优 训练 期 日 
数 来 看 ( 表 3) ,方案 1.1 和 1.3 的 最 优 训练 期 日 数 分 
别 大 于 方案 1.2 和 1.4, 这 是 因为 在 训练 期 日 数 相同 
的 情况 下 ,方案 1.2 和 1.4 的 训练 样本 数 分 别 比方 案 
1.1 和 1.3 多 一 些 , 由 于 有 24 个 预报 时 次 ,前 者 最 多 
可 以 达到 后 者 的 24 倍 ,有 利于 方案 1.2 和 1.4 在 训练 
期 日 数 较 小 时 有 较 多 的 训练 样本 建立 线性 回归 方 
程 。 方 案 13 和 1.4 的 最 优 训练 期 日 数 分 别 小 于 方 
案 1.1 和 1.2, 这 是 因为 方案 1.3 和 1.4 按 照 风 速 是 否 


这 
Mt 
R 
0:00 6:00 12:00 18:00 
时 间 

550 r (c) 方案 1.3 
Z 400F NOO 
w 
R 250 

100 | 


0:00 6:00 12:00 18:00 
时 间 
0.94 — (e) 


0.89 


相关 系数 


0.84 


0.79 


0:00 6:00 
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RI 建 模 方案 | 的 风电 功率 预报 结果 检验 
Tab.3 Test of wind power prediction made 
by the first scheme designs 


统计 量 方案 1.1 方案 1.2 方案 1.3 ”方案 1.4 
最 优 训练 期 日 数 /d 86 64 70 7 
平均 绝 无 误 差 /kW 82 88 48 66 


大 于 $nm:s 对 样本 进行 分 组 ,分 组 后 风电 功率 和 风 
速 之 间 的 相关 关系 更 明显 ,从 而 有 利于 用 较 少 的 训 
练 样本 就 能 建立 拟 合 效果 较 好 的 线性 回归 方程 。 

从 方案 1.1~1.4 的 风电 功率 预报 的 平均 绝对 误 
差 来 看 ( 表 3) ,方案 1.1 和 1.3 的 平均 绝对 误差 分 别 
小 于 方案 1.2 和 1.4, 这 是 因为 一 日 中 风电 功率 与 风 
速 在 09:00 ~ 17:00 时 段 的 相关 系数 小 于 其 它 时 段 ， 
方案 1.1 和 1.3 是 逐 预 报时 次 建立 预报 方程 的 ,在 各 
时 刻 的 预报 误差 都 比较 小 ,而 方案 1.2 和 1.4 是 不 分 
预报 时 次 建立 预报 方程 的 ,在 09:00 ~ 17:00 的 预报 
误差 明显 偏 大 (图 3)。 方 案 1.3 和 1.4 的 平均 绝对 误 


550r b) 方案 1.2 
Z 
et 
R 
00 | | 1 | 1 | 1 J 
0:00 6:00 12:00 18:00 
时 间 
550 r (d) 方案 1.4 
Z 
Mt 
R 
0:00 6:00 12:00 18:00 
时 间 
12:00 18:00 
时 间 


图 3 建 模 方案 I 的 风电 功率 预报 值 与 观测 值 的 逐 时 次 对 比 (a ~ d) 以 及 风电 功率 观测 值 
与 风速 观测 值 的 相关 系数 (e) 


Fig.3 Hourly comparison between observation and prediction of wind power made by the first scheme designs (a ~ d) and correlation 


coefficient between observations of wind power and wind speed (e) 
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一 一 未 按照 风速 量 级 分 组 的 线性 回归 曲线 
一 | 一 按照 风速 量 级 分 组 后 的 线性 回归 曲线 
风速 -风电 功率 散 点 图 


风电 功率 /kW 


UE / mes" 
图 4 风速 一 风电 功率 的 散 点 图 及 线性 回归 曲线 效果 示意 图 


Fig.4 Scatter plot and linear regression curve between wind 


speed and wind power 


差分 别 小 于 方案 1.1 和 1.2, 可 见 按照 风速 等 级 分 组 
后 建立 的 预报 方程 优 于 不 分 组 ,可 用 图 4 来 解释 。 
图 4 任 选 100 次 风电 功率 和 风速 的 观测 资料 作为 样 
本 绘制 风速 -风电 功率 散 点 图 ,并 按照 风速 是 否 大 
F5 ms 将 样本 分 为 2 组 ,分 别 绘制 分 组 前 、 后 风 
速 -风电 功率 的 线性 回归 曲线 ,结果 表明 ,按照 风速 
大 小 对 样本 进行 分 组 之 后 ,风速 一 风电 功率 的 散 点 
更 紧密 地 分 布 在 线性 回归 曲线 两 侧 ,明显 改善 了 风 
速 与 风电 功率 的 回归 关系 。 

图 5 是 对 方案 1.1~ 1.4 的 风电 功率 预报 的 平均 
绝对 误差 的 逐日 对 比 。 从 图 5 来 看 ,方案 1.3 的 平均 
绝对 误差 最 小 ,可 见 , 利 用 一 元 线性 回归 方法 建立 
风电 功率 预报 方程 时 ,引入 风电 功率 与 风速 之 间 相 


300 — (a) 方案 1.1 


200 - 


绝对 误差 /kW 
S 


300 — (9) 方 案 1.3 


200 


绝对 误差 /kW 


关 程 度 的 日 变化 特征 ,以 及 不 同 风 速 量 级 下 风电 功 
率 与 风速 之 间 回 归 关 系 的 差异 性 ,有 利于 减 小 预报 
误差 。 另 外 ,本 文 将 风速 观测 值 假 定 为 预报 值 来 检 
验方 案 1.1~1.4 的 风电 功率 预报 效果 ,不必 考虑 风 
速 预报 质量 随 预 报时 效 的 变化 ,因此 不 必 检 验 较 长 
预报 时 效 的 风电 功率 预报 结果 ,而 取 1~24h 预 报时 
效 是 为 了 考虑 风电 功率 与 风速 之 间 相 关 程度 的 日 
变化 特征 。 

方案 1.3 和 文献 113] 计 算 的 24b 内 的 风电 功率 
预报 的 均 方 根 误差 分 别 在 50 ~ 160 kW 和 200 ~ 400 
kW 之 间 ,前 者 是 在 假定 风速 的 观测 值 为 预报 值 的 基 
础 上 进行 预报 实验 的 ,后 者 则 引入 数值 模式 的 风速 
预报 进行 预报 实验 ,因此 ,不 能 直接 对 二 者 的 预报 结 
果 进 行 对 比 。 要 将 方案 1.3 投 入 业务 应 用 ,还 需要 引 
入 风机 轮 载 高 度 ( 约 80 m) 人 处 的 风速 预报 产品 。80 m 
风速 预报 不 属于 常规 预报 ,例如 ,著名 的 ECMWF 模 
式 不 输出 80 m 风 速 预报 ,但 是 输出 100 m 风速 预报 ， 
可 以 先 将 该 模式 的 100mm 风 速 预报 订正 到 风机 轮 载 
高 度 上 ,再 代入 方案 1.3 进 行 风电 功率 预报 。 
3.2.2 建 模 方 案 开 的 检验 “与 建 模 方案 工 相 比 , 建 
模 方案 荆 的 最 优 训练 期 日 数 和 风电 功率 预报 的 平 
均 绝对 误差 都 明显 偏 大 ,这 与 风电 功率 与 测 风 塔 
风速 的 相关 性 弱 于 与 风机 站 风速 的 相关 性 有 一 定 
关系 。 由 于 建 模 方案 I 和 工 的 研究 资料 不 同 ， 
此 ,不 再 对 建 模 方案 工 与 工 的 风电 功率 预报 结 


800 — (b) 方 案 12 


绝对 误差 /kW 
è 
© 


绝对 误差 /kW 


图 5 方案 1.1~1.4 的 风电 功率 预报 的 平均 绝对 误差 的 逐日 分 布 


Fig.5 Daily variation of mean absolute errors of wind power prediction from scheme 1.1 to 1.4 
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表 4 建 模 方案 [I 的 风电 功率 预报 结果 检验 
Tab.4 Test of wind power prediction made by the 4 结 论 


second scheme designs 


统计 量 方案 21 方案 22 方案 23 GHA 本 文 利用 2016 年 1 月 1 日 ~2018 年 12 月 31 日 
aha ee a e o 陕西 省 某 一 风电 场 区 域内 的 观测 资料 ,分析 了 风电 
平均 绝 无 误 差 /kW 103.58 261.81 261.92 263.2 场 区 域 的 风速 规律 ,并 将 分 析 结 果 融 入 风电 功率 预 
5 报 的 方案 设计 中 ,引进 最 优 训 练 期 方案 ,研究 利用 
线性 回归 方法 建立 风电 功率 预报 模型 的 可 行 性 。 


J 主要 结论 如 下 : 


5j po ee pm A) EPA NRE py Be A 
a es ee eh al rs 现 出 最 大 值 出 现在 夜间 ,最 小 值 出 现在 白天 ,从 低 
倍 WE ’ 了 ,与 仅 以 80 Mm [Ay se 和 风速 作为 预报 因子 层 到 高 层 的 日 变化 趋势 一 致 的 特征 。 该 风电 场 区 


的 建 模 方案 相 比 , 增 加 其 它 高 度 的 风速 作为 预报 因 。 起 各 襄 度 的 风 还 日 杰 化 特征 与 前 人 研究 的 其 它 地 
了 后, 风 电功率 现 报 的 科 沟 绝对 误差 明显 增 大 。 从 区 的 风电 场 有 一 些 差异 ,这 与 各 风电 场 的 下 垫 面条 
二 方案 21~24 的 风电 功率 预报 值 与 观测 值 的 逐日 条 和 低层 大 气 层 结 状态 二 证 益 异 丰 关 


D CE 6) REAN (FINE 5 HEETE, 373R 2.1 KURR O ee eS 
= i ere 日 中 ,二 者 在 09:00 ~17:00 时 段 的 相关 系数 明显 小 
E a 于 其 它 时 段 。 按 照 风速 是 否 大 于 5 ms"' 将 训练 期 
= 文献 13 ERRAMU EAA 2A, T A EE-S UE I 
报 研究 ,发 现 利用 不 同 高 度 风速 作为 预报 因子 比 只 。 的 回归 关系 - 
a eA ee (3) 以 风机 轮 载 高 度 处 的 风速 作为 预报 因子 ， 


要 小 一 些 ,与 本 文 利用 多 元 线性 回归 方法 进行 风电 并 引入 风电 功率 与 风速 之 间 相关 程度 的 日 变化 规 


功率 预报 研究 的 结论 相反 。 这 可 能 因为 风机 功率 。 入 以 及 不 同 风速 量 级 下 风速 与 风电 功率 之 间 回 归 


am ZE Ys eK 42% in| AW — ae ay BRE 5 e : 
PRAT AUE D8 ERS RA — ERRE. 和 最 优 训练 期 短 的 特点 ,可 以 投入 业务 应 用 。 但 是 
1200 厂 (a) 方 案 2.1 1200 三 (b) 方 案 2.2 
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图 6 方案 2.1~2.4 的 风电 功率 预报 值 与 观测 值 的 逐日 比较 


Fig.6 Daily comparison between observation and prediction of wind power from scheme 2.1 to 2.4 
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增加 其 它 高 度 风 速 作为 预报 因子 ,采用 多 元 线性 回 
归 方 法 建立 的 风电 功率 预报 方程 ,预报 误差 较 大 ， 
不 满足 业务 需 3 
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Abstract: Under the background of global warming,disasters caused by extreme temperature changes such as drought, 
heat waves, and cold snaps arebecoming more and more frequent, which poses a serious threat to industrial and agricul- 
tural production, human life and property.Thus, extreme temperature changes have gradually become a widespread con- 
cern. Thus study analyzes temporal and spatial changes of extreme temperatures and explores factors influencing tempera- 
ture indices using daily maximum and minimumtemperature data from 404 stations in northern China for the period 
1960-2017. The study calculated 16 extreme temperature indices using the RclimDex1.0 function of R and investigated 
spatial—temporal variation characteristics of extreme temperature events and influence factors using linear regression, the 
Mann- Kendal test, the sliding t-test,and cumulative anomaly and correlation analysis . The results of the analysisfound 
the following. (1) There was an upward trend in warm extremes and extremal indices and a downward trend in cold ex- 
tremes anddiurnal temperature range; the magnitudes of changes in cold extremes are obviously higher than those of 
warm extremes and the magnitudes of changes in night extremes are higher than those of day extremes.Warming in the 
northwest region is the most obvious, whereaswarming in the northeast region is not obvious. (2) The mutation time of the 
extreme temperature indices mainly occurred in the 1980s and 1990s, and the cold extremes and minimum values of dai- 
ly maximum (minimum) temperature mutated earlier than those of the warm extremes and maximum values of daily maxi- 
mum (minimum) temperature. After mutation,extreme warm events and extreme value events tended to occur frequently, 
whereas extreme cold events gradually decreased in frequency. The extreme temperature indices in the northeast region 
showed the earliest mutation time and the northwest region showed the latest mutation time. (3) Changes in the extreme 
temperature indices also respond to the global warming hiatus of 1988-2012. (4) Most of the extreme temperature indices 
show a close correlation with longitude, latitude, and elevation. (5) The Arctic Oscillation index has the greatest impact 
on extreme temperatures and the most significant effect on the cold index. (6) During 2007-2016, there was a downward 
trend in aerosol optical depth (AOD), and extreme high temperature events increased and extreme low temperature events 
decreased. Most cold indices were negatively correlated with AOD and most warm indexes were positively correlated with 
AOD. This study analyzes extreme temperature variation characteristics in northern China, aims to deeply understand ex- 
treme temperature change trends tohelp carry out meteorological disaster prediction,conduct early warning research,re- 
duce the adverse impacts of extreme temperature events on local production and life, and ensure the safety of agricultural 
production. 


Key words: wind farm; wind power prediction; linear regression method; optimal training period 


